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DEBYE SCREENING IN THE HOT NON-ABELIAN PLASMA 
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The Debye mass is computed nonperturbatively in the deconfined phase of QCD, where 

chromomagnetic confinement is known to be present. The latter defines mo to be mo = co^/aj, 
where cp = 2.06 and cr s = cr B (T) is the spatial string tension. The resulting magnitude of m£)(T) 
and temperature dependence are in good agreement with lattice calculations. 
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1. Introduction 

The screening of electric fields in QCD was 
originally considered in analogy to QED 
plasma, where the Debye screening mass was 
well understood 1 , and the perturbative lead- 
ing order (LO) result for QCD was obtained 



long ago 



,(LO) 



6 



1/2 



gT. An- 



other difficulty of the perturbative approach 
is that the gauge-invariant definition of the 
one-gluon Debye mass is not available. Re- 
cently the calculations of mr}(T) was com- 
puted on the lattice for Nt = 0,2 3 ' 4 ' 5 using 
the free-energy asymptotics. 

The purpose of our talk is to pro- 
vide a gauge-invariant and a nonperturba- 
tive method, which allows to obtain De- 
bye masses in a rather simple analytic cal- 
culational scheme. In what follows we use 
the basically nonperturbative approach of 
Field Correlator Method (FCM) 6J and 
Background Perturbation Theory (BPTh) 
for nonzero T 8 ' 9 > 10 to calculate m£>(T) in 
a series, where the first and dominant term 
is purely nonperturbative 11 

m D (T) = M (T) + BPTh series. (1) 

Here Mq(T) is the gluelump mass due to 
chromomagnetic confinement in 3d, which is 
computed 11 to be M (T) = c D y/<r s (T), with 
ct s (T) being the spatial string tension and 
cd ~ 2 for N c = 3. The latter is simply 
expressed in FCM through chromomagnetic 



correlator 12 ' 13 , and can be found from the 
lattice data 14 ' 15 . Therefore M (T) is pre- 
dicted for all T and can be compared with 
lattice data. 

2. Background Perturbation 

Theory for the thermal Wilson 
loop 

One starts with the correlator of Polyakov 
loops (£(0)L+(r)) 



<L(0)Z+(r)) = 4exp (J^'- T) 



T 



+ - 



iV 2 



■ exp 



F 8 (r,T) 
T 



The first term contributing to the free energy 
F\ of the static QQ-P a i r in the singlet color 
state, the second to the octet free energy F 8 . 

A convenient way to define Debye mass 
was suggested in 3 ' 4 

Fi(r,T) - Fi(oo,T) « -\^l e —o{T)r 

o v 

Accordingly one ends for Fi with the thermal 
Wilson loop of time extension (3 — 1/T and 
space extension r, 

e X p(-pF 1 (r,p)) = (W(r,f3)) 

= i-(trFexp(z«? J ApdzJ). (2) 



i 
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Eq.(2) is the basis of our approach. In what 
follows we shall calculate however not F\, 
which contains all tower of excited states 
over the ground state of heavy quarks QQ, 
but rather the static potential V\ (r, T) , cor- 
responding to this ground state. As a re- 
sult the Debye mass mo is calculated gauge- 
invariantly in terms of the thermal Wilson 
loop. 

Separating, as in BPTh the field into 
NP background and valence gluon field 



A li = B ll 



(3) 



one can assign gauge transformations as fol- 
lows 



B^U+iB^-d^U, 



U+a^U. (4) 



As a next step one inserts (3) into (2) 
and expands in powers of ga^ 7 which gives 
an expansion for the Wilson loop 

(W(r,p)) = (W^(r,/3)) B + 

+ (W^(r,P)) B , a + ..., 
and {W( 2 >} can be written as 



(5) 



j dx fl dy v 



x (trP$(JJ)<a, t ( a; )a„(y)) $(]J))B. (6) 

xy 

Here $(n) an d ^(U) are parallel trans- 
porters along the pieces of the original Wil- 
son loop W(r, /?), which result from the dis- 
section of the Wilson loop at points x and 
y. Thus the Wilson loop (r, (i) is the 
standard loop (r, (3) augmented by the 
adjoint line connecting points x and y. It is 
easy to see using (4), that this construction 
is gauge invariant. 

For OGE propagator one can write 
the path integral Fock-Feynman-Schwinger 
(FFS) representation for nonzero T 

Gp V (x,y) = {a^{x)a v {y)) a 



= l°°ds J(D 4 z): y eM-KW lj (C X y) 
x[P F exp(2ig[ F ap (z(T))dT)) , (7) 

V JO / pv 

<S> ad >(C xy ) = Pcxp(ig f B^Zp) 



where the open contour C xy runs along the 
integration path in (7) from the point x to 
the point y and K = j J" s (i M ) 2 dr. The 



(D z) xy is a path integration measure with 
boundary conditions z m (t = 0) = x^ and 
z M (r = s) = j/ M and with all possible wind- 
ings in the Euclidean temporal direction (this 
is marked by the superscript w). 

We must now average over B^ the geo- 
metrical construction obtained by inserting 
(7) into (6), i.e. 

xy 

m]J)^(C xy M]])) B ee (W xy (r,0)) B . 

xy 

(8) 

For T > T c one can apply the nonabelian 
Stokes theorem, which yields the area law 6 
for distances r >• X g , X g — gluon correlation 
length, X g ~ 0.2 fm 

(W xy (r,f3)) B =eM-Va dj S$), (9) 

where is the area of the space- like projec- 
tion of gluon-dcformed piece of surface S g i, 
and cr^ is the adjoint spatial string tension, 
a adj = C 9 / 4 )^ for SU(3). a s is the funda- 
mental spatial string tension. 

As the result one obtains exactly the 
form containing the gluclump Green's func- 
tion 



B,a 



(*3) 2 C 2 (/) 



x dx 4 dy 4 G 44 (r,t 4 ), (10) 
Jo Jo 

where U = x 4 - y 4 , C 2 (f) = (N 2 C - l)/2N c 
and G 44 (r, t 4 ) is 

G 44 {r,U)= l°° da f(D 4 z): y 
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xexp(-X)exp(-af dj ^). (11) 

Thus the gluon Green's function in the 
confined phase becomes a gluelump Green's 
function, where the adjoint source trajec- 
tory is the projection of the gluon trajec- 
tory on the Wilson loop plane. Gluelump 
is the system when a gluon is connected by 
the string to the Wilson loop plane. String 
is made of chromomagnetic field only. Thus 
the problem reduces to the calculation of the 
gluelump Green's function 

G(x,y) = (x\^2exp(-H n r)\y). 

The lowest eigenvalue of Mo Hamiltonian 
yields the static QQ potential 



V^\r, T) 



N 2 - 1 a s 
2N r r 



-Mnr 



As a result one finds 11 mn = Mq 

(2.82 - 5.08af)yft~ rs 2.06 a 



m D 

This is the central result of our talk. Note, 
that no free parameters are used, a e s s is taken 
from glueball spectrum and BFKL. 

3. Numerical results 

One can now compare our prediction for 
mu(T) = cd \/a s (T) with the latest lat- 
tice data 5 . At temperatures T > Td r ~ 
(1.2 -r 1.5)T C , (where Td r is the temperature 
of dimensionally reduction and physical jus- 
tification for resorting to dimensionally re- 
duced regime at T ~ ^/cr s (T) was given in 
13 ' 16 ) the spatial string tension is chosen in 
the form 14 < 15 

^Jf)=c a g 2 (T)T, (12) 

with the two- loop expression for g 2 (T) 



g- 2 {t) = 26 ln 



where t = T/T c . 

The measured 15 spatial string tension 
in pure glue QCD corresponds to the values 
of c a = 0.566 ± 0.013 and L a = A a /T c = 




Fig. 1. Chromoelectric Debye mass mo/T for 2- 
fiavor QCD (upper lines) and quenched (Nf = 0) 
QCD (lower lines) versus T/T c . Solid lines are cal- 
culated using m D {T) = 2.09 ^/<r s (T), where yJ(T a (T) 
corresponds to Eq.(12) with Nt = 2 for upper solid 
line and Nt = for lower solid line. Dashed lines are 
calculated using Eq.(14), Nt = 2-upper line, Nt = 
0— lower line. The lattice data are from O.Kaczmarek 
and F.Zantow, Phys. Rev. D71 114510 (2005). 

0.104 ± 0.009. On the Fig. 1 arc shown lat- 
tice data 3 and theoretical curves for the De- 
bye mass in quenched (Nf — 0) and 2-flavor 
QCD. Solid lines correspond to our theoret- 
ical prediction, mr>(T) = cd \/ 's(r), with 
co = 2.09 and for a/ct s (T) we exploit the 
parameters c a = 0.564, L a = 0.104. The 
upper solid line is for the Debye mass in 2- 
flavor QCD, and the lower - for quenched 
QCD. We note that in computing moiT) us- 
ing (12), (13) all dependence on Nf enters 
only through the Gell-Mann-Low coefficients 
bo and b\ . For comparison we display in the 
Fig. 1 dashed lines for m£>(T)/T, calculated 
with a perturbative-like ansatz 5 . 



Latt ( 
l D 



N 



(T)=A Nf ^l+'-fg(T)T. (14) 

Comparison of our theoretical prediction 
with the perturbative-like ansatz shows that 
both agree reasonably with lattice data in 
the temperature interval T c < T < 5T C ; the 
agreement is slightly better for our results. 
At the same time, in (14) a fitting constant 
is used Apf, ~ 1.5, which is necessary even 
at T/T c ~ 5. We expect the accuracy of the 
first approximation of our approach is around 
10%, taking also into account the bias in the 
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definition of for the gluelump. The tem- 
perature region near T c needs additional care 
because i) the behaviour of \J u s (T) deviates 
from the lattice data and ii) contribution of 
chromoelectric fields above T c which was ne- 
glected above. 

4. Conclusions 

We have studied Debye screening in the hot 
nonabelian theory. For that purpose the 
gauge-invariant definition of the free energy 
of the static QQ-pair in the singlet color state 
was given in terms of the thermal Wilson 
loop. Due to the chromomagnetic confine- 
ment persisting at all temperatures T, the 
hot QCD is essentially nonperturbativc. To 
account for this fact in a gauge- invariant way 
the BPTh was developed for the thermal 
Wilson loop using path-integral FFS formal- 
ism. As a result one obtains from the ther- 
mal Wilson loop the screened Coulomb po- 
tential with the screening mass correspond- 
ing to the lowest gluelump mass. Apply- 
ing the Hamiltonian formalism to the BPTh 
Green's functions with the einbein technique 
the gluelump mass spectrum was obtained. 
As a result, we have derived the leading term 
of the BPTh for the Debye mass which is the 
purely nonperturbative, 

m D (T) = c D Vct s (T), c D ^2. 
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